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1.0 INTRODUCTION 


The Interplanetary Program to Optimize Simulated Trajectories (IPOST) is 
intended to support many analysis phases, from early interplanetary feasibility 
studies through spacecraft development and operations. The IPOST output 
provides information for sizing and understanding mission impacts related to 
propulsion, guidance, communications, sensor/actuators, payload, and other 
dynamic and geometric environments. 

Much of the overall architecture for IPOST has been derived from the Program to 
Optimize Simulated Trajectories (POST). Indeed certain POST parameters and 
capabilities have been incorporated into IPOST to aid in POST-IPOST user 
compatibility. IPOST has extended trajectory capabilities to target planets and 
other celestial bodies with intermediate and velocity correction maneuvers. 
IPOST capabilities and limitations are summarized in Table 1>1. 


FEATURE 

CAPABILITY 

Optimization method 

Explicit (Master/subproblems), Implicit (collocation) 

Optimization algorithm 

NPSOL 

Optimization 

parameter 1 * 1 

AV magnitude, mass, time, . . . 

Maximum controls 

25 (Master), 45 (subproblems), 1700 (collocation) 

Control parameters* 

Values of event criteria, AV, arrival conditions, thrust, 

Maximum targets 

25 (Master), 45 (subproblems), 1700 (collocation) 

Target parameters* 

Time, position, velocity, orbital conditions, . . . 

Targeting method 

NPSOL, Newton-Raphson, special Onestep 

Sensitivity matrix 

Finite differencing, analytic for special interplanetary 
targeting 

Maximum events 

100 

Event criteria* 

Time, distance, speed, closest approach, . . . 

Event activities 

Info, impulsive AV, launch, orbit insertion, mass 
jettison 

Maximum maneuvers/ 
subproblems 

15 


Conic, Onestep, Multiconic, Encke, Cowell, implicit 

Planetary bodies 

Sun, nine planets, Earth's moon, any user-defined 
bodies 

Ephemeris 

Analytic, precision (JPL) 

Trajectory 

perturbations 

Central body, perturbing bodies, radiation pressure, 
J2, aerodynamics, thrust 

Input/Output frames 

Ecliptic or planet equator, Mean 1950 or Mean 2000 

| * User selectable | 


Table 1 - 1. IPOST Features/Capabilities 
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IPOST, along with members of its family, such as POST and IPREP, can analyze 
and support almost every activity associated with space exploration (Figure 1- 1). 
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Figure I - 1. Space Exploration Infrastructure 
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IPOST is event driven. That is, the user defines a sequence of events which are 
executed in the simulation process. The events can be triggered by different 
criteria, such as absolute or relative time, distance from a body, or propellant 
consumption. At the event times, various activities can be initiated or 
terminated, such as employing a different thrust steering law, changing 
trajectory propagators or propagation step size, performing an impulsive delta 
velocity maneuver or jettisoning a probe or stage. 

The time period between two contiguous events is called a phase. Trajectory 
propagation takes place in each phase. Five types of propagators are available 
(listed in order of increasing accuracy and decreasing computational speed): 
Conic, Onestep, Multiconic, Encke, Cowell. Propagator selection depends upon 
user needs, such as simple fast simulations for parametric feasibility analysis, or 
precision detailed trajectories to support subsystem design. 

IPOST can run a single trajectory simulation or it can run multiple simulations. 
For multiple simulations, one can run a parametric scan and/or an optimization 
mode. The search mode will vary one parameter, such as planetary arrival time, 
over a specified interval and increment size, and perform a simulation (or 
optimization) for each search parameter value. 

The optimization mode will optimize a user cost/objective function, such as 
mBTimiim mass that can be placed in a desired orbit, subject to user-specified 
constraints. The constraint variables, such as periapsis altitude or orbital 
inclination, are called dependent variables or target parameters. The parameters 
which are free to vary, such as maneuver delta velocity (AV), are called 
independent variables or control parameters. As part of, or instead of, 
optimization, trajectory targeting can be performed. In this case, there is no cost 
function and the IPOST problem reduces to finding a set of control parameter 
values that meet specified target parameter conditions. 

Generalized targeting and optimization uses the Stanford NPSOL algorithm. For 
certain types of problems, a trajectory decomposition method is available. There 
is a master optimization process which requires that the trajectory be divided into 
legs or sub-problems. Each sub-problem is an optimization problem in itself, 
containing controls, constraints and an (optional) objective function. A special 
application of decomposition is the Interplanetary Targeting and Optimization 
Option (ITOO). This technique uses ansdytical partials generated during nominal 
trajectory propagation to determine minimum AV (or mass) trajectories, usually 
for gravity assist (swingby) missions. 

In addition to the classic method of explicit optimization, there exists an option to 
perform implicit optimization using the collocation method. In this case, each 
phase is divided into independent segments which are allowed to vary subject to 
intersegment continuity and the equations of motion. Optimization using 
collocation is less sensitive to faulty initial guesses, but requires much greater CP 
time than explicit optimization to achieve the same level of accuracy. 
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IPOST input is via three namelists: $TOP, $TRAJ and $TAB. $TOP contains a 
description of the targeting and optimization problem. It must be input first. 
$TRAJ contains data that describes each mission event/phase. It must follow 
$TOP , and there must be one $TRAJ for each event. $TAB is used to input 
tabular data such as thrust vs. time or drag coefficient vs. mach number and 
angle of attack. Input and output units are metric. 
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2-0 APPTJCATTONS 

Examples of mission applications are shown in Table 2-1. They illustrate some 
of the diverse IPOST model capabilities, including lunar, interplanetary, orbital, 
thrusting (impulsive, low, high), and gravity assist. 

The typical IPOST application is usually in the form of run sequences which 
compare various mission options. These analysis threads build on each other, 
culminating in a reference mission which is used to support detailed system 
design and analysis. Table 2 - 2 illustrates representative mission threads which 
exercise key IPOST capabilities. 

For example, in the Comet Rendezvous thread, each case is actually several runs 
of IPOST to generate parametric data such that mission decisions and 
refinements can be made. The sequence of multiple runs per case feeds each 
succeeding case with each case becoming more realistic in terms of model 
fidelity, and encompassing more system objectives and constraints. 

The first step in the Comet Rendezvous thread is to define mission requirements 
(necessary conditions) and goals (desired conditions), as well as any known 
constraints, such as technology status. In this thread, a mission requirement 
would be to successfully rendezvous with a specific comet in a specific time frame. 
A mission goal might be to collect data on planets or bodies that are encountered 
during the interplanetary trajectory from Earth to the comet. A mission 
constraint might be the availability of a Cesium ion thrust engine powered by 
solar arrays, which provide limits on power/thrust levels and on specific impulse. 

The first case is an approximate impulsive AV solution using a Venus gravity 
assist. This examines coarse energy requirements, benefits of gravity assist, and 
optimum mission opportunities (launch-arrival dates with payload/launch mass 
sensitivities). The mission may be analyzed as separate phases, e.g., Earth to 
Venus, Venus to comet approach, rendezvous and stationkeeping. 

The second case models a low thrust mission using a single thrust segment with 
variable steering and variable throttle. Implicit optimization (collocation) is used. 
This recognizes the coarseness of the initial guess, and provides rapid solution 
searches. A determination is made whether available technology is sufficient to 
provide the required payload mass at comet rendezvous. 

As mission knowledge evolves, the third case introduces multiple coast/thrust 
segments. These added degrees of freedom provide more flexible, and more 
realistic, mission solutions. The optimization method can be implicit or explicit 
(with master-subproblems), depending upon how many, and what level, of 
mission decisions need to be made. This would include interplanetary and close 
encounter geometries, flight times, subsystem performance, etc. 
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TEST MISSION PLANETARY 

CASE DURATION BODIES ACTIVITIES SIMULATION 
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The final case provides an end-to-end precision optimized reference trajectory for 
system analyses and subsystem design support. Science and mission objectives 
can be assessed with a high degree of confidence. 


• Lunar Mission Thread 

LI Earth departure to lunar orbit with patched conic 

L2 Space Station to moon with free return 
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• Voyager II Thread 
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• Comet Rendezvous with Solar Electric Propulsion Thread 

Cl Approximate impulsive AV (DSM) with Venus gravity assist. 
C2 Single thrust segment, simple collocation 
C3 Multiple thrust segments, complex collocation 
C4 Explicit optimization comparison 

• Human Mission to Mars Thread 

Ml EM launch/arrival date search, simple aerocapture 
M2 MVE return leg optimization, Earth orbit capture 
M3 Mars surface ascent to orbit (MAV design) 
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3.0 SAMPLE CASES 

The following sample cases illustrate detailed mission applications of IPOST. 
These are by no means intended to cover all capabilities nor are they realistic in 
every detail, but they do provide meaningful examples for constructing and 
understanding mission applications. There are 8 case discussions, four include 
both IPOST input and output, and the last four cases include only IPOST input. 

The first sample case is a Voyager 2 mission which illustrates a Master- 
Subproblem formulation. The second sample case is a manned Mare mission 
that shows a complete round-trip. The third sample case is a Lunar mission 
which departs from Space Station Freedom and terminates with an insertion into 
lunar orbit. The fourth sample case illustrates optimization using collocation for 
an Earth-Jupiter mission. 

The fifth case illustrates trajectory simulation with no targeting or optimization. 
The sixth case is a version of Voyager 2 with subproblem optimization. The 
seventh case shows a classic Hohmann transfer problem using finite thrust and 
collocation. The eighth case is a low thrust mission to Jupiter. 
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3.1 VOYAGER 2 

The Voyager 2 case illustrates a master subproblem formulation and multiple 
planetary encounters. Only the Earth-Jupiter-Satum phase of Voyager is 
performed. The simulation has 8 events, starting in Earth park orbit and ending 
with a Saturn flyby. Total impulsive delta V is minimized in the master problem. 
The two subproblems target each of two legs. Earth to Jupiter and Jupiter to 
Saturn. 

In setting up an IPOST problem, the trajectory simulation is defined first, as 
opposed to the optimization process, because it describes the primary mission. 

The first $TRAJ namelist is event 5. This utilizes parameters such as S/C mass 
and propulsion characteristics. The initial date of July 31, 1977 precedes the 
actual Voyager launch date. The S/C is placed in a circular orbit about Earth, 
and trajectory propagation will use conic, or two-body, equations of motion. 

The second event (#10) is triggered by a flight time of 20 days. The 1STEP 
propagator is activated in connection with activation of the LAUNCH mode. 
Hence, Voyager orbits the Earth for 20 days and then is impulsively injected onto 
an escape hyperbola. Using flight time as a control parameter would allow 
variation of initial launch date. For 1STEP, the primary body is the Sun and 
Earth is the secondary body. 

The third event (#15) is initiated alter a flight time of 20 days from the launch 
event. At this point, the S/C is well outside of the Earth's sphere of influence. The 
secondary body is now defined as Jupiter for 1STEP propagation. The reference 
body for trajectory calculations is also set to Jupiter. 

At event 20, the triggering criteria is mission time. The intended value is Jupiter 
closest approach time which is not specified explicitly, but indirectly through the 
optimization process. A conic propagator is used with Jupiter as the primary 
body. 

One day later, at event 23, an impulsive trajectory correction maneuver is 
executed. 1STEP is reactivated as the propagator as the vehicle flies away from 
Jupiter. The combination of gravity assist and midcourse correction will set up 
the trajectory for the flight towards Saturn. 

The fifth event (#25) occurs 300 days after the midcourse maneuver. The 
secondary body for 1STEP is switched to Saturn, as is the reference body. A 
general rule of thumb is that planetary sphere's of influence for 1STEP (and 
Multiconic) are about 20 days for small inner planets (Mercury to Mars) and 200 
days for large outer planets (Jupiter to Neptune). 

The next event (#30) is Saturn's closest approach, on August 25, 1981. It occurs 
1467 days after Earth park orbit escape. 
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A final event (#90) is used as a "dummy" event, and coincides with the previous 
event. This event is needed because final event computations are done only on the 
"minus" side of the event. The namelist specification of "NONE" means no more 
input follows; that is, this is the end of IPOST problem specification. 

We now return to the optimization definition, or $TOP namelist. The NPSOL 
algorithm is used with finite difference perturbations being calculated internally 
by NPSOL. This is the normal mode, as opposed to the user specifying 
perturbations. A maximum of 10 iterations is allowed for master problem 
optimization and 600 iterations for each subproblem targeting process. A 
subproblem targeting iteration is typically much shorter than a master problem 
iteration because the latter must solve all subproblems successfully. 

The master problem controls are Jupiter closest approach time (defined as the 
criterion of event 20), B dot T and B dot R of the Jupiter approach, or incoming, 
hyperbola (defined as the 2nd and 3rd dependent variables of subproblem 1). In 
addition to the initial guesses for control values, the upper and lower bounds, and 
weighting values are important inputs. These define the performance manifold 
and often mean the difference between problem convergence and divergence. The 
objective function is the sum of all delta V magnitudes, which in this case 
corresponds to the launch/escape maneuver at Earth plus the Jupiter midcourse 
correction. 

The two subproblems are defined next. The first subproblem ends at event 20 
(Jupiter closest approach) and the second subproblem ends at event 30 (Saturn 
closest approach). IPOST automatically assumes that subproblems are non- 
overlapping. Both subproblems use a Newton-Raphson targeting technique, as 
opposed to subproblem optimization with NPSOL. 

The control variables for subproblem 1 are the V-infinity vector of the Earth 
departure hyperbola. Controls for subproblem 2 are the AV components after 
Jupiter flyby. As in the master problem, important inputs are the control initial 
guesses (USUB), bounds (INDSLB and INDSUB), and weightings (WGTS). 

The constraint or target variables for subproblem 1 are time from periapsis, B dot 
T and B dot R at Jupiter. For subproblem 2 the constraints are time from 
periapsis, B dot T and closest approach distance at Saturn. The use of B dot T as a 
control parameter (with loose bounds) at Saturn affects orbit inclination as well as 
which side of Saturn the S/C flies by. B dot T and B dot R are often used because of 
their stability in the targeting and optimization process. The choice of constraint 
parameters for subproblem 2 reflect termination of the mission at Saturn flyby. 
For the actual Voyager mission, which continues on to Uranus and Neptune, the 
constraint parameter types at Saturn would have been identical to those at 
Jupiter. 

Because IPRINT was not input, the default value of IPRINT = 0 is used. This will 
result in only summary information of the master problem, plus the final 
trajectory, being output. Except for well tested production runs, it is 
recommended that more detailed levels of print be exercised. 
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The first page of output, after the Namelist echo, summarizes IPOST input, 
including initial conditions, event summary, targeting and optimization 
definition, master-subproblem structure, and NPSOL options. 

The next page of output completes NPSOL parameter definition and then prints 
an iteration summary for each major (master level) iteration. This includes the 
objective function value (sum of delta V) and other conditions of the optimization 
process. Of some interest are the number of objective function evaluations 
(NFXJN) the condition number of the Jacobian (COND T). This problem does 
not have any nonlinear constraints, only control bounds, but if it did, the value of 
the nonlinear constraint norm is displayed, as well as the constraint values. 

As the iterations progress, certain key parameters should be monitored. These 
are the objective function value (which should decrease), the condition numbers of 
the Hessian and Jacobian (which should remain small), and the value of the 
nonlinear constraint norm (which should decrease). Also of interest are the 
convergence indicators at the right side of the summary. When all these flags are 
"T", then successful optimization has been achieved (according to NPSOL). One 
cautionary note is that reliance on a few indicators, such as the convergence 
flags, can be misleading. It is important to examine all measures, including the 
final solution. 

The iteration summary continues until page 4 of the output, where a maximum 
iteration limit has been reached (INFORM = 4). NPSOL is exited, and the final 
trajectory is displayed. Conditions at each event are output, including the minus 
and plus sides of each impulsive maneuver. At the end of each subproblem, a 
final iteration summary is displayed. In this case, Subproblem 1 convergence is 
summarized between trajectory blocks for event 15 and event 20 (Jupiter dosest 
approach), and Subproblem 2 is summarized just before event 30 (Saturn closest 
approach). 

The last solution has a total AV of 7.27 km/s, reduced from the initial solution of 
7.67 km/s. Jupiter dosest approach time is 711.9 days from launch, and B dot T, B 
dot R were sdected as 1898747 km, 130453 km. 

Finally, a master problem summary output is given which includes the number 
of iterations and CPU time. A valid solution may or may not meet all convergence 
criteria. Only the user, with adequate engineering experience, can make that 
judgment. 
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indslb = -20., -20., -20 . , -10 . , -10 . , -10 . , 

indsub = 20., 20., 20., 10., 10., 10., 

wvus = 20., 20., 20., 10., 10., 10., 

pertsb = .00001, .00001, .00001, 3*. 00001, 

wgts = .01, .180849d7, ,129268d6, 15., .01, ,206130d6, ,2d6 



depsvl - -.0001, .180849d7, . 129268d6, 0001 206130d6, 200000 . dO 
depslb * -.0001, . 180849d7, . 129268d6, 0001 2061 30d6, 1 99999. dO 
depsub * -.0001, .10O849d7, . 129268d6, -. 0001 ,. 2061 30d6, 200001 . dO 
depstl « .00005, .1, .1, .00005, 100., 100.' 
wvsnlc - 1., . 180849d7, . 129268d6, 1 . , . 206130d6, .2d6 
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Variables 

7. 118782E+01 0 1 . 049873E+00 0 1.014660E+00 
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Nonlinear objective value = 7.270559E+00 


tables 

7.119208E+01 0 1.049746E+00 0 1.007911E+00 
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Nonlinear objective value * 7.269251E+00 

Values of the constraints and their predicted status 



2.2E-01 38 7*269194E+00 003 7.7E+00 
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state relative primary body: sun 

x =- 0 . 974184351 E +09 y = 0 . 229413788 E +09 z = 0 . 328818838 E +08 
vx =- 0 . 126002267 E +02 vy ’- 0 . 125491850 E +02 vz = 0 . 770292069 E +00 
sma — 0 . 256300375 E +10 eccen «= 0 . 1293484 10 E +01 inc = 0 . 25698994 3 E +01 
argp = 0 . 351055672 E +03 anlong *= 0 . 1 19694447 E +03 meaan = 0 . 71 9129706 E +01 
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argp = 0. 130483133E+03 anlong = 0 . 15461 7110E+02 meaan = 0. 216645061E+03 
initial conditions for phase 90 
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esn = 90.000 fesn= 90.000 

t ime= 2 . 4 4484250D+06 
normal termination 
cpu = 43.180 seconds 
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3.2 MANNED MARS 


The manned Mars case is launched from Earth park orbit, performs a Venus 
gravity assist, orbits into a high elliptic Mars orbit, remains in orbit for two 
months, then launches back to Earth where the S/C inserts into an elliptical 
holding orbit. The simulation has 13 events with a total mission time of 610 days. 
This mission is a pure targeting case with no master controls, constraints, or 
objective function. All targeting is done on the subproblem level. Each 
subproblem is full rank (number of controls = number of constraints). 

To understand this case we begin with the simulation, or first $TRAJ input. At 
event 5, the initial conditions are specified, such as a S/C mass of 1000000 kg and 
an initial circular park orbit of 6878 radius. 

The next event (#10) performs an impulsive launch maneuver, .5 days later, to 
insert the S/C on an Earth escape hyperbola. The 1STEP propagator is activated to 
account for 3 body motion (Earth-Sun-S/C). The S/C will leave Earth and begin its 
journey toward Mars, via Venus. 

At event 20, 40 days after park orbit escape, the secondary body for 1STEP is 
changed to Venus (from Earth). The reference body (IDBODY) for target 
calculations is also switched to Venus. 

Events 25 and 30 are simultaneous, but contiguous events. Event 25 is triggered by 
a mission flight time of 190 days. This is intended to correspond to Venus closest 
approach. Event 30 performs an impulsive midcourse correction. 

At event 35, 40 days after Venus closest approach, the secondary body and 
reference body are changed to Mars. 

Event 40 and 45 are simultaneous contiguous events. Event 40 corresponds to 
Mars closest approach, 350 days after Earth escape. Event 45 performs an 
impulsive orbit insertion into a highly elliptical Mars orbit whose inclination is 
the same as the approach hyperbola. The propagator is switched to a simple 
conic. At this point subsequent events could have simulated a Mars surface 
descent and subsequent ascent. However, for test case simplicity, only the 60 day 
orbital phase is simulated. 

At event 50, after a two month stay in orbit, the S/C is inserted into an escape 
trajectory towards Earth. The 1STEP propagator is reactivated with Mars as the 
secondary body (and the Sim as primary). After 40 days, event 60 switches the 
secondary body and target reference body to Earth. 
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The final 3 events (#65, #70, #75) are all simultaneous contiguous events. Event 
65 corresponds to Earth closest approach after a total mission time of 610 days. 
Event 70 performs an impulsive orbit insertion into a highly elliptical holding 
orbit about Earth. The orbit might be typical of a staging orbit which transfers 
crew and material to and from the low Earth Space Station orbit. Event 75 exists 
to allow conditions after the orbit insertion (before event 75) to be used as output 
values. 

We now return to the optimization formulation, or $TOP. No optimization 
algorithm was selected and only subproblem targeting will be done. Forward 
differences are specified for the finite differencing construction of the Jacobian. 
The maximum subproblem targeting iterations is 300. 

The 3 subproblems all use the N ewton-Raphson targeting algorithm. The 
subproblems correspond to Earth-Venus, Venus-Mars, Mars-Earth legs. The 
first subproblem varies escape V-infinity to achieve desired B dot T, B dot R, TFP 
at Venus closest approach. The second subproblem varies Venus midcourse delta 
V to achieve desired closest approach distance, time, and B dot T at Mars. The 
third subproblem varies Mars escape V-infinity to achieve desired closest 
approach distance, time, and B dot T at Earth return. 

The output for this case is straightforward. 

The final trajectory is displayed, corresponding to each event in the specified 
simulation input, along with the subproblem summary output (achieved by 
setting ISUB (i) = 1 in $TOP). Of particular note are the conditions before and 
after each maneuver. The final S/C mass at Earth return is 70404 kg, reduced 
from an initial mass of 1 million kg in Earth park orbit. 
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Interplanetary post simulation, version 2.18 , dated 03-05-90. problem no. 
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final conditions for phase 35 
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initial conditions for phase 50 after launch maneuver 
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vx - 0.190043391E+02 vy — 0.818599076E+01 vz - 0.610586295E+00 speed - 0.207014049E+02 fpa — 0. 666317000E+02 
sma “”0 . 930133391E+03 eccen - 0.488496640E+05 Inc - 0.174136625E+02 argp - 0.843450076E+02 anlong - 0. 162094779E+03 
meaan — 0.647756922E+07 truan — 0.666327767E+02 tfp — 0.587928657E+02 rperi - 0 . 454357735E+08 vperl = 0 .207016606E+02 
Incoming Asymptote 
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total cpu time - 55.390 seconds 


3.3 LUNAR ORBITER 


This case examines flight from a low Earth orbit (LEO) into a low lunar orbit. It 
has a master-subproblem formulation, and minimizes total AV. 

The simulation starts with a specified Cartesian state vector at Earth (which 
corresponds to circular orbit of 6878 km radius). An impulsive maneuver is 
identified. The conic propagator is activated with Earth as the primary body. At a 
mission time specified as an optimization control, an impulsive midcourse 
correction is applied, at which time the primary and reference bodies are 
switched to the Moon. 3 days after Earth orbit escape, the S/C arrives at lunar 
closest approach and inserts into lunar orbit. 

The optimization defines NPSOL as the optimization algorithm using finite 
difference perturbations computed internally by NPSOL. 50 iterations are allowed 
for the master problem and 200 iterations each for the 2 subproblems. The master 
control parameters are the Earth escape delta V components and the mission 
times of the midcourse maneuver and lunar closest approach. 

The first subproblem varies midcourse AV to achieve desired hyperbolic 
conditions (closest approach distance, inclination, time from periapsis) at lunar 
arrival. The second subproblem varies lunar insertion AV, applied at hyperbolic 
periapsis, to achieve desired lunar orbit conditions (periapsis and apoapsis radii, 
and inclination). 

The output echoes the namelist $TOP and all $TRAJ namelists. A summary of 
the problem is given, and execution output begins with a summary of each 
iteration of each subproblem for master iteration 0. Following the first converged 
iteration (13) of subproblem 2 is the summary of master iteration 0. The sequence 
repeats for all master iterations. Only master iterations 0, 1, 2, and the final 
iteration (13) are shown. After the iteration summary is displayed, the final 
trajectory simulation is output. 

The total AV was 3.968 km/s, reduced from 4.058 km/s for the initial guess. Note 
that the iteration summaries display the weighted objective function value. The 
midcourse maneuver time has been moved from 2.39 days to 1.50 days, the lower 
control bound. Lunar dosest approach time has been moved from 3.00 days to 3.5 
days, the upper control bound. 
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interplanetary post simulation, version 2.18 , dated 03-05-90. problem no 
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3.4 EARTH TO JUPITER (COLLOCATION) 


The following sample case is a portion of the Voyager mission, Earth to Jupiter 
phase, to illustrate collocation. Of particular interest is the difference between 
explicit optimization such as the test case described in Section 3.1, and implicit 
optimization, or collocation. 

The mission starts (event 10) with a S/C mass of 2213689 kg and an input state 
corresponding to a post park orbit escape bum. Ideally, the initial state should 
produce a trajectory which leaves Earth and flies by Jupiter with the desired 
encounter conditions. This would be true for an accurate explicit propagator. 
With collocation, trajectories are defined as segments between events, and by 
Hermite cubic polynomials between segment end points, or nodes. In this case, 
the number of segments per phase is set to 1 for all phases. This represents a 
rather poor trajectory approximation, unless the events are placed very close 
together. More will be said about this later. An impulsive AV maneuver is placed 
at this first event (ideally the AV should be zero), and at every event except the 
final event. 

i 

The next event (#11) occurs 14 minutes later (.01 days), where another impulsive 
AV is performed (again, ideally equal to zero). 

There are 6 more events in the case (#12, 13, 17, 18, 19, 20) which are planned at 
mission times of 0.21 days, 2.21, 666.1, 684.1, 685.6, and 686.1 days. The last event 
(#20) ideally corresponds to Jupiter closest approach. 

The placing of event times is an attempt to normalize accelerations within each 
phase. If there is a large difference in accelerations, or trajectory curvature, 
within a segment (= phase for this case), then the Hermite polynomials are very 
poor representations of a real trajectory. Inaccurate Hermites will cause 
divergence in the optimization process. 

In this case, the optimization specifies "colloc" with "autoperts". The MXITOP 
array indicates that Jacobian rescaling occurs immediately after it is initially 
calculated by finite differencing. A maximum of 150 iterations is allowed after 
rescaling. 

The mission control variables are the AV components of the impulsive maneuver 
at each of the first seven events. Other input control variables for the collocation 
formulation are specified by the U(22) array elements. These correspond to the 
Cartesian state vectors at each event from #11 through #20, a total of 7 
(components per state) x 7 (events) = 49 elements. These state vectors were 
obtained from a previous explicit optimization run. In practice, the state vectors 
could have been obtained by a "WAG” simulation run since collocation can be 
robust. 


75 


The constraints are TFP, B dot T and B dot R at Jupiter closest approach (event 
20). Note that lower and upper bounds for the constraints may contain more than 
3 elements each. This is not a problem because the collocation process overrides 
all non-mission constraints, that is, specifies upper and lower bounds for 
constraints introduced by the collocation process. The collocation process also 
computes weightings for the non-mission constraints. 

The test case output first echoes namelist inputs. Then a summary of NPSOL 
parameters is displayed. The problem begins with computation of the Jacobian. 
IPOST will only allow non-zero elements of the Jacobian to be numerically 
computed. The explicit zeroing out of null elements saves a considerable amount 
of CP time because a collocation Jacobian will be a sparse matrix. In this case, 
there are a total of 126 controls and 108 constraints. Of the 13608 Jacobian 
elements, about 9600 are theoretically non-zero. 

The first (before rescaling) iteration (0) summary follows, Of critical interest is 
the Jacobian condition value (COND T). A value of 3.E7 before rescaling and 3000 
after rescaling indicates a poorly conditioned Jacobian. The independent (control) 
and dependent (constraint) parameter values are displayed after each iteration 
summary. 

The next iteration takes place and the resultant independent and dependent 
parameter values are displayed after the iteration 1 summary. Noteworthy items 
indude the Jacobian condition of 40000 (a significant worsening), norm of the 
nonlinear constraints of 8.84 (also degraded), and an objective (total AV) of 4.40 
(internal units). 

Intermediate iterations are omitted to save space. The final converged iteration 
(#140) is displayed. The Jacobian condition number has improved slightly, but is 
still bad. The norm of constraint violation is very small, indicating constraints 
have been met. However, this has occurred at the expense of an objective value 
increase to 5.19. 

The final solution is displayed after convergence. In this case, a whopping 51.9 
km/6 total AV is required. This large AV indicates an aberrant solution. Figure 3 
- la is an x-y plot of the final trajectory. This points out an interesting 
phenomenon. Even though NPSO L has found an optimal, and a "feasible", 
solution, as indicated by die TTT convergence flags, the solution is not realistic 
from an engineering perspective. A local optimum has been found. Sometimes 
the mathematical robustness of collocation is deceiving. 

The problem in this case is the number of segments per phase. When they are 
increased to two segments per phase, then a solution illustrated in Figure 3 - lb 
occurs. The total AV in that case is 14 km/s. The input file shown can be used for 
two (or more) segments per phase by changing NSEGPH to more than 1. With 
COINTP = 'LINEAR', IPOST will interpolate linearly between the event states to 
supply the internal node states. 
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Unfortunately, collocation is not as simple as always increasing the 
number of segments per phase. As with any optimization problem, 
particularly one with high dimensionality such as collocation, the 
local performance manifold is strongly determined by weightings 
(independent and dependent parameters), bounds, and initial 
conditions. These can be very tricky with the large dimensions 
associated with collocation, and may produce misleading solutions. 
Caveat emptor! 
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X coantifiaac (kn) 


a) 1 Segment/Phase, AV = 52km/s 



•7—C* -W+A4 W*« *4#+0i -2e*04 -W*« 0 U«« 2^M IMI 

X coordinate (km) 


b) 2 Segments/Phase, AV a 14 km/s 
Figure 3-1. Ecliptic View of Earth to Jupiter Mission 
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indxi (4 ) = 4,5,6, 

indvr (4 ) = * dvx' , ' dvy' , ' dvz' 

indph (4) * 3*11, 

u<4) = 3*0., 

wvu (4) = 3*0.01, 
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critr - 6htimrfl, 

value « 0 . 66613102e+03 
lpbody (1) « 5, 

idbody » 5, 

mantyp *= 'impuls'. 
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parameter octal decimal 

event criteria data - 531b 345 
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//NPFEAS// The maximum violation is 6.08E+00 in constraint 157 
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/ /NPFEAS/ / The maximum violation is 3.54E-08 in constraint 196 
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altp =-0. 30668594 3E+0 4 vinfm = 0.426531710E+02 vinfx » 0 .306846113E+02 vinfy *=-0 .296154901E+02 vinfz * 0 . 818761245E+00 
bmag = 0 . 352357178E+04 btheta = 0 . 120650090E+03 bdt =-0 . 179629473E+04 bdr = 0 . 303131710E+04 hypta = 0 . 8644 19184E+02 
c3 = 0 . 181 929299E+04 ra =-0 . 424 4524 06E+02 dec = 0 . 109990522E+01 altit « 0 . 814 436176E+07 
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date = 8 22 1977 5.04 julian = 2443377.71000000 tdurp = 2.00000000 critr = timrfl propid = lstep 

prlmid = earth secid = earth idbody = earth frame * ecliptic epoch = mean 2000 

timrfl = 0 . 22100000 0 E +01 
state relative to idbody: earth 

x = 0 . 58649157 7 E +07 y =- 0 . 565807241 E +07 z = 0 . 153429146 E + 06 radius = 0 . 815073990 E +07 scmass = 0 . 593693257 E +01 
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2£ OTHER CASgS 


The following cases contain only the IPOST input, and related discussion. The 
input corresponds to executable IPOST runs, and are shown to expand the 
envelope of IPOST usage. 


3.5.1 TRAJECTORY SIMULATION 

The input file described here performs explicit trajectory propagation. A Cowell, 
or numerical integration of the equations of motion, is emphasized. There is no 
targeting or optimization. This rather simple case can be used to compare 
propagators or effects of forces, to perform sensitivity analyses, to generate initial 
guesses, such as for collocation rims, or to debug a simulation. 
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3,5,2 CO M PLE TE YQ YAC rE Rn 


This case is an extension of the Voyager 2 case (Section 3.1). The mission 
continues beyond Saturn, and includes flybys of Uranus and Neptune. In 
addition, the Newton-Raphson full rank subproblem targeting is replaced by 
NPSOL optimization for each subproblem. The master problem remains the 
same, although central differencing with input perturbation values is used to 
form the Jacobian. 

In the subproblem setup MODELT is set to 'SUBOPT' for each subproblem. Each 
subproblem still has 3 controls and 3 constraints, but now target error is 
minimized except for the first subproblem where AV is minimized. 
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indxsi = 1,1,1, 2,2,2, 

indsvr = 5hvinf x, 5hvinf y, Shvinf z , 3hdvx, 3hdvy, 3hdvz 

indsph = 3*10, 3*23, 

usub = 2.38521,9.44533,3.19406, 3*0., 

pertsb = 3*1 .d-7, 3*1 .d-5. 



indxsi (7) = 3,3,3, 4,4,4, 

indsvr(7) = 3hdvx, 3hdvy, 3hdvz, 3hdvx, 3hdvy, 3hdvz 
indsph { 7 ) = 3*33, 3*43, 
usub(7) = 3*0., 3*0., 
pertsb (7) = 3*1. d-6, 3*1. d-5. 



in 





T5 










o 





m 





00 





r- 





H 





i 





VO 





T3 

m 




r* 

'P 




cm 

vo 




VO 

m 




ON 

ov 




m 

vo 




* m 

vo 




u 

T5 





x *h 

f-H 




X o 

«. o 




m o 

-w o 




«. o 

AJ 

Ll O 

ft) 




T3 1 

X 

a i 

Ll 




4= 

X 




m vo ■ 

in m 




- *o - 

* T> * 




a on «h 

CL • 




VM o • 

lW ON O 




AJ VO % 

aj in fH 




X h h 

X m » 




m cm o 

n r- rH 




«-• o 

CM O 




% ■ o 

* • O 




M «. • 

Ll *0 




T5 r- 

TJ vo O 




J3 *0 * 

X T3 • 




x m •— i 

X vo 




*m vo o 

m vo *> 




oi % m • 

* in • 




% jj %m « 

% aj *-h 




eg TJ o h 

vOon •. 




.jjmno 

» J2 IT) H • 


k H 


CM X * • • 

X * ■ «-H 


«. m 



n<M . . 


*P - 


% * t— 1 *— 1 

«. *i H rl 


c r- 


a *o o 

m cl * o o 


0) - 


%VM o o o 

vVM o o o 


<-h r- 


HiJCMOO 

cn aj o o 


qj r- 


*X * • O 

*X * • o 

* 

U ON 


*-h m m i 

n n n i • 

in 

* »— ♦ 


ll H ll II H 

ll ll ll ll ll 

ii 

II II 



^ 





r- r-* r- r- r- 



10 





AJ 

TJ M£HH 

"O 


X 

0 

to > a > aJ 

to > a > jj 

aJ 

u 

O' 

x to W 10 W 

x to to io in 

t“N C 

o a) 

L| 

ti a cl a cl 

'O CL CL CL CL 

qj 0 ) 

CL AJ 

qj 

C t) t) 0) 1) 

C H) HI Cl ID 

u > 

0) (4 

AJ 

•H XI T3 T3 T3 

*H TJ *P *P *P 

A» 4I-H*0 



</> 

<A 


c 

c 

0 

U 

o a 



o 

o 


u 

•H % 

a U 

— < C VOiO * 

O O T3 T3 O 

0) * - O • • 00 

.> ro m h(n v 

II II II H II II II 


E >1 >1 to aj 

qi T3 TJ a M 10 

M O O O « 3 

**-t X 43 u E l» -h 
•o a^j a u x a 

•H *H *H *H « 4J 10 


o 


c m 

o VO 
U to 
% VO 

II II 


3 

a 

c 



i-H 





VM 


X 



Li 

CL 

u 



E 

a) 

c 

• . . 


*H o. 

AJ » 


m m lo 


«. AJ • 

to m 

q} 

vo vo • 


0X0 

«— i *. 

«— i « 

m m co 

ft) 

a 

qj 

HVO (M 

«. o 

■. CM 

VO VO CM 

II II II 

n ll 

II II 

II II II 

O 





to 





ft) 





X 



CL 

a 

AJ 

AJ Ll (U 

CLT3 

>iX 

■H qj 

f'l Li 

C U 3 

O O 

X U 

Li O 

qj qi 

(0 -rl H 

Li X 

C C 

o a u 

M 0) 

> Li qj 

a cl 

qj 

a qj c 

AJ 

ft) U > 

-rH -rH 

E *H 

Lt Ll -H 


a u u 


</> </> 

a u o 


122 


change idbody 
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critr = Shtdurp, 

value = 100., 

mantyp = 1 lmpuls' 

dvx * 3*0. dO, 



critr - 6htimrfl 

value = 3099., 
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3.5.3 HOHMANN TRANSFER 


The Hohmann transfer problem is a classic orbital mechanics solution: The 
optimal two impulse solution for transferring from one circular orbit to another 
circular orbit in the same plane is to perform the impulses 180° apart in the 
direction of travel (for raising an orbit). 

This case uses finite bums and collocation. Although the case assumes a low 
enough thrust such that "impulse" is not a good approximation, and only one 
segment per phase is assumed, the Hohmann transfer solution in principle is 
validated. 

The first bum occurs in the phase between the first event (# 10) and the second 
event (# 20), followed by a coasting phase, and then the second bum starting at 
event 30 and ending at the last event (# 40). Simple thrust tables are input, in this 
case constant thrust (5000 newtons). 

There are 7 mission controls: Thrust direction relative to velocity, throttle, and 
stop time of the first bum, and start time, thrust direction relative to velocity, 
throttle, and stop time of the second bum. Invalid guesses of the states are input 
at each event. There are a total of 49 additional, collocation introduced controls. 
The constraints are Cartesian position at the end of the second bum. Final mass 
is maximized. 
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hohmann transfer orbit test case 
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initial conditions for phase 30 
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3.5,4 LOW THRUST TO JUPITER 

This case illustrates a nuclear powered low thrust mission from Earth to Jupiter. 
There is one long continuous burn. The control parameters are thrust throttle 
level, and pitch and yaw thrust directions (UVW relative). Constraints are TFP, 
and closest approach distance at Jupiter arrival. S/C mass at arrival is 
maximized. 
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Galileo 1989 VEEGA trajectory 
im * 'npsol'. 
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